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Primordial black holes with a narrow mass range are regarded as a nonrelativistic fluid component 
with an equation of state for dust. The impact of the black hole evaporation on the dynamics of 
the early universe is studied by resorting to a two-fluid model. We find periods of intense radiation 
reheating in the initial and final stages of the evaporation. 



I. INTRODUCTION 

The interest in primordial black holes (PBHs) dates 
back to Zel'dovich and Novikov |Q, Hawking g|, and 
Carr and Hawking ||, who concluded they should be 
produced by density fluctuations in the primordial cos- 
mic fluid, leading to a wide mass spectrum for these ob- 
jects. The earlier their formation, the lower their mass. 
Since then there has been a manifold of proposals for 
the production of PBHs: Inhomogencities triggered in 
the inflationary period of cosmic expansion Q, first or- 
der phase transitions , bubble wall collisions (t| , the 
Gross-Perry- Yaffe mechanism of quantum gravitational 
tunneling from hot radiation j|] (see also [g| and [[To||), 
quantum pair creation (based on the nonboundary pro- 
posal of the wave function of the universe) during infla- 
tion 1 1 1 1 , quantum fluctuations in hybrid inflation ]T^ ] 
(see also @), and the decay of cosmic loops [|4), to 
mention just a few. 

This interest in PBHs is understandable as they may 
have had a notable impact on several areas of current in- 
terest, such as baryogenesis [ p~5| , dark matter frl6| ], galac- 
tic nucleation [p"7| , and the reionization of the universe 
p8[ . In addition they can constrain the spectral index of 
primordial density perturbations fl9| , |2C| . 

A natural outcome of the copious production of black 
holes in some of the mentioned scenarios is that, sooner 
or later, a point will be reached where the energy density 
of the PBH component will significantly contribute to 
the energy density of the cosmic medium. For example, 
phase transitions based on grand unified theories or 
bubble wall collisions in extended inflation may lead to 
a PBH dominated universe. Also, in scenarios of hybrid 
inflation there may exist a regime in which the black 
holes will dominate the total energy density soon after 
the end of inflation [Jl2| . This situation may last for some 
time before the contribution of the latter component to 



in some respect plays the role of an equilibrium 



the total energy density becomes negligible because of its 
evaporation via Hawking's process. (This stage has been 
termed "the binary phase" and studied in detail j2lj for 
PBHs produced by the Gross-Perry- Yaffe mechanism.) 

PBH formation processes from phase transitions, e.g. 
in connection with extended inflation, are characterized 
by a very narrow mass spectrum (see, e.g., |fl2[ and refer- 
ences therein). In such a situation the assumption of only 
one mass m {BH) for all members of the population may 
be considered a reasonable approximation. On the other 
hand, a black hole of mass m (BH) may thermodynami- 
cally be characterized by a temperature T {BH) tx fn~ BH) 
p3[ . Consequently, with a single mass population of 
PBHs one may associate a single temperature as well. 
This suggests a picture of the PBH component in which 

T <.. 

fluid temperature 

On the basis of this interpretation it is the target of 
the present paper to study the cosmological dynamics 
during PBH evaporation within a model of two interact- 
ing and reacting fluids. The evaporation process in this 
model is pictured as a decay of the PBH "fluid" into a 
conventional fluid. We discuss the reheating of the latter 
as a consequence of the evaporation and show that one 
may describe this process in terms of a bulk pressure of 
the cosmic medium as a whole. The general formalism is 
then applied to the dynamics of a Friedmann-Lemaitrc- 
Robertson- Walker (FLRW) universeand the back reac- 
tion of the PBH evaporation on the cosmic scale factor 
is investigated. 

To establish a transparent picture we introduce the 
following simplifying assumptions from the outset: (i) the 
PBHs can be thought of as particles of a perfect fluid with 
the equation of state for dust, (ii) both fluids share the 
same four-velocity, (iii) all the PBHs start evaporating 
at the same time. Once these assumptions are accepted 
one can apply the formalism of interacting and reacting 
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fluids to the mixture of both components (say radiation 
and the PBH fluid) |J]. Notice that if the PBHs did 
not evaporate, the evolution of the scale factor and the 
energy densities of both components would directly be 
given by the well-known results of Jacobs for a mixture of 
relativistic particles and dust [g5| . However the fact that 
the PBHs radiate makes a big difference with respect to 
the rather traditional case of a mixture of radiation and 
dust with conserved fluid particle numbers. 

A closely related work of Barrow et al. |^(| considers 
the evaporation of PBHs with an initial power-law mass 
spectrum. The black holes in that paper do not consti- 
tute, however, a fluid in the sense discussed here (such 
kind of fluid picture looks less natural for a nonsingular 
distribution of masses). In general, the dynamical be- 
havior of both configurations is different, especially with 
respect to the existence of solutions for which the energy 
density in PBHs is in equilibrium with radiation. 

This paper is organized as follows. Section II recalls 
the general hydrodynamical formalism for two interact- 
ing and reacting perfect fluids |24|j2^ ] which previously 
has been applied to the reheating phase at the end of the 
inflationary expansion p8| , p9| |. Section III specifies this 
framework to a mixture of thermal radiation plus pri- 
mordial black holes and derives a general formula for the 
corresponding entropy production density, which is then 
used to discuss specific features of the evaporation pro- 
cess, especially a "reheating" of the radiation. The evolu- 
tion equations for the energy densities of each fluid cou- 
pled to the Friedmann equation are solved numerically. 
Section IV introduces an effective one-component model 
of the cosmic medium and describes the PBH evapora- 
tion process in terms of an effective bulk viscous pressure 
of the system as a whole. We present a general formula 
for the latter, which is evaluated explicitly in section V 
for the initial stages of the evaporation process. Finally, 
section VI summarizes the results of this work. Units 
have been chosen so that c=kB = h = G = l. 



II. THE GENERAL TWO-FLUID MODEL 

The energy-momentum tensor T %k of the cosmic 
medium is assumed to split into two perfect fluid parts, 
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where is the particle number density. We are inter- 
ested in situations where neither the particle numbers nor 
the energy-momenta of the components are separately 
conserved, i.e., conversion of particles and exchange of 
energy and momentum between the components are ad- 
mitted. The balance laws for the particle numbers are 
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n (A)T(A) 



(4) 



where O = u\ is the fluid expansion and T (A) is the rate 
of change of the number of particles of species A. There 
is particle production for r.^. > and particle decay for 
T (A) < 0, respectively. For T (A) = we have separate 
particle number conservation ("detailed balance"). 

Interactions between the fluid components amount to 
the mutual exchange of energy and momentum. Conse- 
quently, there will be no local energy-momentum conser- 
vation for the subsystems separately. Only the energy- 
momentum tensor of the system as a whole is conserved. 

Denoting the loss and source terms in the separate bal- 



ances by t l , we write 
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implying 
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P(A) + © (p W +P (A) ) = U at a (A) , 



\i = -hff 



(A) 



(5) 



(6) 



(7) 



All the considerations to follow will be independent of 
the specific structure of the t l . In other words, there 

(A) 

are no limitations on the strength or the structure of the 
interaction. 

Each component is governed by a separate Gibbs equa- 
tion: 



T <A) ds (A) =d ^+P<A) d : 
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with 



where T {A) and s {A) are the temperature and entropy per 
particle of species A, respectively. Using the balances (^) 
and (^), one finds for the time behavior of the entropy 
per particle 



(A) =P(A)UU + P(A) h 
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where p {A) denotes the energy density and p {A) the equi- 
librium pressure of species- A particles. For simplicity 
we assume both components to share the same four- 
velocity u l , normalized so that u l Ui = —1. The tensor 
h lk = g tk + u l u k projects any vectorial quantity on the 
hypersurface orthogonal to u l . The particle flow vector 
N l of species A is defined as 

(A) L 



AjV) =«aC ~ (P {A) +P i A ) ) T 
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With nonvanishing source terms in the balances for n (A) 
and p {A) , the change in the entropy per particle is differ- 
ent from zero in general. 

The equations of state of the fluid components are as- 
sumed to have the general form 



V) \ n (A)' T (A)/ 



P { A) 0V)> T (A)) > ( 10 ) 
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i.e., particle number densities n (A) and temperatures T (A) 
are regarded as the basic thermodynamical variables. 
The temperatures of the fluids are different in general. 

Differentiating the last of the relations (|ni|), using the 
balances ([|) and (|^) as well as the general relation 



d P(A) _ P(A) + P(A) T (A) d P 
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that follows from the requirement that the entropy is a 
state function, we find the following expression for the 
temperature behavior ppj ] : 

t = -T (Q-T ) -^L 

(A) (A) V w x (A)) a 

U V(A) 



+ 



(A) ' P{A) 



dp ( A } /dT w 
Here we have used the abbreviations 

9p (A) _ dp w /dT w dp w _ ( dp w 
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There is no corresponding condition, however, for the 
particle number balance as a whole. Defining the integral 
particle number density nasn = n (1) + n (2) , we have 



n + 9n = nT = n (1) T {1) + n (2) T ( 



(17) 



r is the rate by which the total particle number n 
changes. We do not require T to be zero since total par- 
ticle number conservation is only a very special case. 



The entropy per particle i 
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where p [A) is the chemical potential of species A. Intro- 
ducing the expression (O) into Eq. (If4|) yields 



u a t a 

n r + w 

(A)\a rp (A) (A) ' rp 



S a — ^ {A) 
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For the total entropy production density S a a = . a 
S a .„ we obtain 



Both the source terms T {A} and u a t® in Eqs. (Q) and 

(^J) back react on the temperature. The numerator of the 
second term on the right hand side of Eq. (|fj) coincides 
with the right-hand side of Eq. (Q), i.e., the correspond- 
ing terms disappear in the special case s. A . = 0. 

For r, .. — u„t a — and with = 3a la, where 

(A) a (A) ' 1 

a is a characteristic length scale (which in the homoge- 
neous and isotropic case coincides with the scale factor 
of the Robertson- Walker metric) , the equations of state 
P<2) = n (2)^(2): P(2) — 3" (2 )T'(2) l° r radiation (subindex 
(2)) reproduce the well-known XL, ~ a -1 behavior. With 
p (1 , = n (1) T m , p w = n m m+ §n fl) X (1) one recovers 
XL. ~ a -2 for matter (subindex (1)). 



The entropy flow vector S a is defined by 



S W ~ n (A) s (A) u ° 
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and the contribution of component A to the entropy pro- 
duction density becomes 



° (A y,a - (A) (A) (A) 
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where relation (^) has been used. 

According to the balances (||) the condition of energy- 
momentum conservation for the system as a whole, 
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The condition S a a = requires the well-known equilib- 
rium conditions (see, e.g., p2| (chapter 5)) 



P(l) P {2 ) I "V) ^(2) I 

as well as r = 0. 



(21) 



III. BASIC RELATIONS FOR AN 
EVAPORATING BLACK HOLE COMPONENT 

In this section we apply the formalism outlined above 
to a mixture of massless radiation (component 2) and 
evaporating PBHs (component 1). We assume all the 
black holes of the same mass m {BH) , and that they evap- 
orate solely into radiation particles, though this is admit- 
tedly a rough approximation |33|. The PBH component 
is treated as a nonrelativistic fluid with 



P. 



and the energy density 



= 



(22) 
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(BH) 



n (BH) m (BH) 



The number N^ BH) of PBHs in a comoving volume a 3 , 
= n (BH) a 3 1 is not preserved and, according to Eq. 



N, 



(|4|), we may write down a balance equation for the cor- 
responding PBH number flow vector N* — n [BH) u l , 



N 



y-)TI = ?? I 

^'"(BH) (BH) (BH) 



(24) 



^From Eq. (|23j) the black hole energy balance becomes 
[cf. Eq. (§)] 



Op,„„. = u a t a 

" (BH) (BH) 



(25) 



with 



u a t a 

a (BH) 



P(BH) 

The entropy per black hole is 

S,„„, = 47TTO 2 



(BH) 



m , 



(26) 



(27) 



while the black hole temperature is related to its mass 
by the well-known formula |23| 
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(BH) 



8ttiti, 
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The temperature T. BH . is attributed to each PBH indi- 
vidually, i.e., it is not a conventional fluid temperature. 
With this basic conceptional difference in mind it is nev- 
ertheless possible to associate standard fluid-type quan- 
tities to the single-mass PBH component. Specifying the 
general relation ([t8]) to the case at hand, a chemical po- 
tential may be associated with the black hole component 
by 



T s 

(BH) (BH) 



1 

2 m (BH) 



(29) 



Together with relation (|28|) the expression ( |23| ) for the 
PBH energy density fits into the general structure given 
in Eqs. (fffl. The corresponding partial derivatives are 



d P(BH) 

dT„ 



(BH) 



d P(BH) 
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Furthermore, the relationship 



d ^(BH) 



1 m 



(BH) 



2 T„ 
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(31) 



holds. It is essential that the general temperature law 
(12) is applicable as well and specifies [cf. Eqs. ( |22| ) and 
(26)] to 



u a t a 

" (SB) 
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(BH) " (BH) 



d P(BH)l dT (BH) 



-T, 



(BH) m. 



(32) 



(23) which is consistent with formula (|2£|). The BH temper- 
ature behavior T {BH) oc m ~ BH) may be regarded as a spe- 
cial case of the general fluid temperature law (TH\). This 
provides the main motivation for our fluid approach to 
PBHs. On the other hand, this approach relies on the 
possibility of treating the black holes as an ensemble of 
nonintcracting particles. The validity of the latter as- 
sumption can be shown by a simple Newtonian argument 
pl| . The relative velocity of neighboring black holes due 
to the Hubble expansion is 

Vexp = dt (v 1 ") 3 ) • 

With the help of Eq. (E4T) we obtain 
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d» t „ = — n "'" I 1 '— - 

cxp q (BH) \ Q 



r, 



3 (BH) 

for this quantity. Since for a spatially flat FLRW universe 



e /8tt , 



3 V 3 
is valid, we find 



v exp 



\bS\/y ( n (BH) m (BH) +P(2)) (l - 
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The velocity required for a BH to escape the (Newtonian) 
gravitational pull of its nearest neighbor is 



2m,„„,n, nlf = n 1 ^ 3 1 /2n 



(BH) '"(BH) 



(BH) ""(BH) 



Comparing the expressions for v exp and v esc , we find that 
v exp > v esc , provided 



\ 




> 1 



The latter relation is always satisfied for I\ SH . < 0, the 
case of main interest here, but also for a certain range 
of values r,„„. > 0. We conclude that the behavior of 

(BH) 

the individual BHs is primarily determined by the expan- 
sion of the universe and not by their mutual interaction. 
This justifies their treatment as a nonintcracting particle 
species. 

Let us assume from now on that the source terms on 
the right-hand side of Eq. (||) cancel among themselves, 
i.e., that the entropy per particle of component 2 is pre- 
served. This isentropy condition amounts to the assump- 
tion that the fluid-2 particles at any stage are amenable 
to a perfect fluid description. With s (2) = 0, via Eq. (||) 
equivalent to 

U ^2) = (P(2) + P(2)) r (2) > ( 33 ) 

the expression ( |l2| ) for the temperature behavior simpli- 
-,ly: 



fies considerably: 
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(- 



[ (2)) d P(2) 



Combining Eqs. (|16|) and (|33|), one has 

u a t a = ~u a t a = - (p, oS + p, n .) T, 

a (SB) a (2) V (2) ^(2) J { 

and we find the expression 

P(BH ) 



P(2 



(34) 



(35) 



(36) 



for the rate of change of the fiuid-2 particle number. Neg- 
ative rates r <Bff) and rh, Ba , /tti. bh , (evaporation) imply 
a positive T (2) and vice versa. 

The contribution of the black hole component to the 
entropy production density is [cf. Eqs. (|l3) and (|24|)] 



na 
(BH) -a 



with 



according to Eq. (£ 
tributes with 



Qa 



(37) 



= 2^^ , (38) 

m (BH) 

T). The fluid- 2 component con- 
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(39) 



at t = in. To obtain the second relation ([y]) we have used 
Eqs. ( p7|) and (f28|). Obviously, this case requires r < 
in order to satisfy S? a > 0, i.e., decay (T, < 0) of 
the PBH component. In other words, the formation of 
PBHs (T (BH) > 0) is thermodynamically forbidden if the 
temperature of the PBHs coincides with the temperature 
of the ambient radiation, i.e., for T {BH) — T {2) . 

Let us now consider the second term in Eq. ( |40|) . This 
contribution, which vanishes for T {BH) (to) — T {2) (to), 
represents the entropy production density in case the 
number of PBHs is constant, i.e., for F {BH) — 0. Restrict- 
ing ourselves momentarily to a fixed BH number, there 
are two possibilities for the subsequent evolution of the 
system. For T < T immediately after t = t a pos- 
itive entropy production requires r n {BH) /m {BHj < 0, i.e., 
PBH evaporation. This process will be studied in some 
detail below. For T,„„. < XL, on the other hand, the 
condition S a n > is only satisfied for to, D „, /m. „ „, > 0, 

,Cl " (B H) i {B H) ' 

i.e., for mass accretion. Becoming larger, the BHs also 
become colder [cf. Eq. (|28|)] and their growth seems to 
continue indefinitely, i.e., until they have eaten up the 
entire universe. However, such a kind of scenario is ther- 
modynamically impossible as can be seen by the follow- 
ing argument. The essential point is that the BH mass 
accretion back reacts on the temperature of the ambient 
radiation. For an accretion rate m,„,„ lm,„„. > the 

. . ( B H ) I {BH ) 

rate ( |36| ) is negative (recall that r (BJf) = in the present 
discussion). The cooling rate (|34|) of the radiation tem- 
perature then becomes 



Assuming now the second component to be radiation 
with p {2) = p {v /'S and [i (2) = 0, the overall entropy pro- 
duction density (|20|) becomes 



P{BH)^ {BH) 



2T„ 



{BH) 



T T 

^(BH) -*-(2) 



(40) 



where we have used relations ( |39| ) and ( |36| ) . This formula 
for the entropy production density is still completely gen- 
eral and, within the fluid picture, holds for PBH for- 
mation (r ( > 0) and mass accretion {m (BH) > 0) 
as well as for the evaporation process (T (BH) 



rn, 



< 0). 



{BH) 

< and 

For T,„, > 2T, n „. the second law favors 

{2) {BH) 



the formation of PBHs (T {BH) > 0) and mass accretion 
(m (BH) > 0), while for T < T (BH) the reverse processes 
dominate. 

Let us assume that at some initial time to the black 
hole temperature coincides with the radiation tempera- 
ture, i.e., T (t Q ) = T, 2) (to)- Under this condition the 
second term in Eq. (ETj) vanishes and the entropy pro- 
duction density reduces to 



S" 



1 P{BH) p 

2 J 1 (BH) 



-n s r 

{BH) (BH) (BH) 



(41) 



T 



I [ I 3 P (BH) m {BH) 



4 P(2 



It is obvious that from some time on the temperature 
T (2) will cool off faster than T (BH) [cf. Eq. @]. Con- 
sequently, T (2) will approach T {BH) . As soon as T (2) 
has fallen below T {BH) , mass accretion stops since for 
T (2) < T (BH) the rate r n {BH) /m {BH) has to be negative in 
order to guarantee a positive entropy production, i.e., the 
process now proceeds in the reverse direction. We con- 
clude that the second law of thermodynamics forbids a 
catastrophic growth of the PBHs. These considerations 
may also be regarded as a justification of the "initial" 
condition T [BH) (to) — T (2) (to) for the evaporation pro- 
cess. 

Now we continue discussing the case with Tbh dif- 
ferent from zero. The situation S a a — generally cor- 
responds to a "detailed balance" at t — to, i.e. both 
reactions, the formation of PBHs and their evaporation, 
proceed at the same rate such that there is no change in 
the net numbers N, nrrt and N,„ . This is equivalent to 

(BH) (2) 

r (B H, (*o) = m (B H) (*o) /m (BH) (to) = T (2) (to) = 0. 

For t < to the formation of PBHs and accretion, i.e., 
r (Bff) > 0, rh (BH) /m (BH) > and T (2) < are thermody- 
namically favored. For t > to the PBH evaporation with 
r (B ff) < 0' 1 ^ l {BH)/ m {BH) < an d L (2) > dominates. 
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We are interested here in the second part of this pro- 
cess, i.e., in the evolution of the universe for t > to- 
We will not discuss issues of PBH formation (see the in- 
troduction for possible mechanisms) but simply assume 
that at some early time a considerable amount of cosmic 
matter was in the form of PBHs of the same mass and 
that this component subsequently (e.g., until the epoch 
of nuleosynthesis) decayed. We will investigate the ther- 
modynamic aspects of this evaporation process and its 
implications for the cosmological dynamics. 

According to Eq. ( |36| ) negative values of the rates 
r (BH) and m (BH) /m (BH) imply a positive T (2) . The pro- 
duction rate T (2) may either be larger or smaller than 
the expansion rate 6. For r, 2) < 6 the fluid temper- 
ature decreases according to Eq. (pi[), while the BH 
temperature increases according to Eq. (|32|). It follows 
that T [BH) > TL, at t > to. The evaporation process 
will continue since T {2) < T {gH) requires T (BH) < and 



(BH) 



l m {BH) < to guarantee S? a > in Eq. 



For r. 2 . > 0, however, hypothetically realized e.g. by 
a large initial ratio p (BH) /p (2) , the fluid temperature in- 
creases. If this increase is smaller than the increase in 
T (BH) we have again T ( < T (B and the PBH evapora- 
tion goes on since it remains thermodynamically favor- 
able (S.° > 0). But an increase in T stronger than that 
in T (BH) results in a fluid temperature which is higher 
than T For T (2) > 2T (BH , 



duction ( |40| ) requires r > and fh, BH) /m (BH) > 



a positive entropy pro- 
0, 



implying a quick transition to a negative r <2) , i.e., the 
process confines itself. A strong "reheating" of the fluid 
will stop the evaporation and reverse the process. Now, 
the second law requires PBHs to be formed out of the ra- 
diation and accrete mass. A negative T (2) , on the other 
hand, will make T (2) subsequently decrease [cf. Eq. (53)] . 
If T (2) has fallen below T {BH) , the evaporation process 
may continue (see the discussion following Eq. (El|)), 

We conclude that the PBH evaporation is a self- 
confining process and, consequently, the accompanying 
"reheating" of the radiation is limited on general ther- 
modynamical grounds. 

To study the corresponding dynamics of the two- 
component system in more detail we resort to the decay 
law Wi 



A 



m , 



(42) 



where A may be taken as a positive-definite constant to 
be determined by the condition m {BH) (to + r) = 0, with 
r the lifetime of a Schwarzschild black hole that obeys a 
law of the type ([|2|) from its formation up to its final dis- 
sapearance, leaving no remnant behind. We assume that 
all PBHs start evaporating at the same time and that 
accretion of the ambience matter on the black holes can 
be ignored (|,|2l[] . In reality Eq. ( f42| ) is an approximation 
to Hawking's law since A is not rigorously constant over 
the entire radiation period, but depends on the number 
of particle species emitted by the black hole, whence it 



slowl y i ncreases with the inverse of the black hole mass 
(see jgH and J34|). For the rate th, Bli) /m. BH) one finds 



(BH) 



1 ( m (BH) ^O) 



m (BH) 3r V m (BH) W / 

Integration of this equation yields 

t-to 



1H) (*) = m ,BH) (*<>) 



1 - 



1/3 



(43) 



(44) 



The last expression allows one to write the evaporation 
rate (EH) as 



1 



3rl 



(45) 



The number decay rate of the PBHs is determined by 
their inverse lifetime, i.e., 



(BH) 



(46) 



Formulas (Ej^) and (^) hold if the PBHs are purely de- 
caying with the reverse process entirely suppressed. 

In order to describe the process of formation and sub- 
sequent evaporation of PBHs in full detail, one needs ex- 
plicit expressions for the rates r (BJJ) and w (BH) /m (BH) 
through the transition period from predominant PBH 
formation to the evaporation phase, where r (BJJ) and 



Jm. 



switch from positive to negative values. 



(bh) and r (B ff)' respectively, are 



While our general framework up to Eq. ( }40|) is able to 
cover this most general case, the simple expressions ( E5| ) 
and ( [f6"l ) for m {BH) /m 
valid only for pure evaporation. 

To simplify our description we will assume the inter- 
val between the time at which V IBm changes its sign and 
the time at which the expressions (45) and (E^) are valid 



to be negligibly small, so that it may be justified to ap- 
proximately identify both times. Under this condition 
the simple expressions ( ff5| ) and ([h|) may be used for the 
whole range t > to- As a consequence, we have to deal 
with a nonzero decay rate already at t = to and the "ini- 
tial" entropy production rate is positive, according to Eq. 

©• _ _ 

Applying the simple rates ( |45| ) and (Eq) from t Q on will 

provide us with a transparent picture of the evaporation 
process, although it implies the assumption of an equal- 
ity of T {BH) and T (2) at a time to somewhat later than 
the "real" beginning (T (BH) — rh, 
evaporation process. 

Assuming p ( = p (2) /3 in Eq. 
t is 



IH )/ m (BH) = 0) ° f tllC 

(p6|), the production 



rate F {2) at t 



r w (to) = £ 



P(BH) ^O) 
P(2) (*0) 



(47) 



where (3 is the initial ratio of the PBH energy density 
to the radiation energy density. Given the PBH lifetime 



G 



t, the initial production rate is determined by the ratio 
of the energy densities. According to Eq. ( [34}) this rate 
fixes the subsequent behaviour of T (2) . It is obvious from 
Eq. (H), that for a rate T (2) = 3H - 3m (BH) /m {BH) , 
where H = ® is the Hubble parameter, the fluid tem- 
perature TL, increases at the same rate as T {BH) [cf. Eq. 
([52|)]. Since evaporation is only possible for T {2) < T {BH) , 
this is the highest rate for which there is evaporation. 
A still higher rate interrupts the evaporation and ther- 
modynamically favors a temporary PBH formation phase 
until T (2) w T (BH) again (see the discussion following Eq. 
([fl|)). Combining this maximum rate at t — to with the 
rate Jl5| ) at t = to, provides us with the condition for 
evaporation, 



(3 < 3H t + 1 



(48) 



where Hq is the Hubble parameter at to- This condi- 
tion limits the initial ratio (3. For r < Hq 1 the PBHs 
evaporate within one Hubble time and (3 is restricted to 
(3 < 4. If the BH lifetime is much larger than the Hubble 
time, i.e., r 3> H^ 1 , the initial PBH abundance may be 
larger, i.e. p {BH) > p (2) corresponding to T (2) > T (BH) , 
is possible at t = t . 

Notice that these are restrictions following from the 
thermodynamics of the decay process. We did not take 
into account here any limits on the fractional abundance 
from PBH formation mechanisms which might be more 
restrictive than the present ones. 

reh 

We may define a "reheating" temperature T as the 
maximum temperature for the radiation component by 
T™ — 0. According to Eq. ( |34| ) the corresponding con- 
dition is T (2) = 9. With p {2) = 3n (2) T (2) and Eq. @, 
i.e., assuming the reheating to proceed in the initial phase 
of the evaporation process, we find 



1 U (BH) m (BH) 

3 n {2) 9r 



(49) 



The energy balance ( p5| ) with the "sink" (£6|), applied 
to the evaporation process, becomes 



e P( 

bh) - Pi BH) 



111, 



(50) 



with the rates (^) and (^). For the radiation compo- 
nent we have 



P(2) + O P(2) = -P(BH) 



(BH) 



(BH) 



(51) 



With respect to "equilibrium" solutions found by Bar- 
row et al. |2(|, in which the ratio of the energy den- 
sity in the PBH component to the energy density of 
the ambient radiation remains constant, we address the 
question whether a corresponding behavior is possible in 
the present context. To this purpose we have to study 



whether a relation p {BH) = (3p m with (3 
patible with the equations (pOh and ($1 

P(bh) — 0P(2) m either Eq. (|5Q ) or Eq 
with a relation 



1 (3 



4/3 + 1 



-e 



const is corn- 
Introducing 
provides us 



(52) 



with r, 



--J[ T (BH) +rh (BH) /m (BH} ] [cf. Eq. ©]. 
It follows that such kind of equilibrium solutions re- 
quires the rate T (2) (and, equivalently, the sum r. BH . + 
m (ifl) l m ( B H)) t° be proportional to the expansion rate 
O. But this is clearly incompatible with the expressions 
([45| ) and (^). We conclude that equilibrium solutions 
do not exist in our two-fluid model. This is at variance 
with Barrow et al. (see equation (24) of Ref. J2(|) who 
(for a different PBH configuration) found a correspond- 
ing equilibrium behavior for the evolution of P {BH) and 

P(2)- 

Integration of equations (|5(]) and ( |5l| ) with the rates 
(E5|) and (B^) provides us with the expressions 



P(BH) (t) - P(BH) ^O) ^ 



t-t 



1/3 



exp 



t-t 



(53) 



for the energy density of the black hole component, and 



P W tt) = P(2) (t ) + 



«0 , 1 P(BH) ^O) a 



t 

t 4 _ t-t 

3 



cr 

dt a (t) 







[1 



t-tr 



1 2/3 



exp 



t-t 



(54) 



for the radiation energy density, respectively. The dy- 
namics of the entire system is determined by the energy 
densities ([53]) and ( |54] ) together with the field equations 
for a (t) . In the case of a spatially flat FLRW universe 
the latter reduce to the Friedmann equation 



3 72 = 8lT (P(BH) +P(2)) 



(55) 



We have numeric ally integrated the integro-diffcrcntial 
set of equations (jHq), (j|J) and (^5|) and shown the re- 
sults for p {BH) and p {2) in figures 1 and 2, respectively. 
We have included the possibility of no black hole at all 
((3 = 0) for the sake of completeness and comparison 
with the evolution of a uncontaminated FLRW radiation 
universe. The impact of the evaporation process on the 
expansion rate may be understood in terms of an effec- 
tive bulk viscous pressure of the cosmic fluid as a whole 
as will be discussed in the following sections. It is well 
known that a bulk pressure tends to increase the expan- 
sion rate p5| , ^6[ . 

Denoting the number of relativistic particles by N {2) — 
n (2) a 3 , the following general relations hold for the radia- 
tion component through the whole process: 



7 



oc 



N i/3 

(2) 



T (2) « 



ATl/3 

(2) 



.1/3 



(56) 



implying p {2) 

Toward the end of the evaporation, i.e. for t — > to + r, 
tha rates Th, BH) /m, BH) and r <2) diverge, leading to ex- 
plosive particle emission in the final stages of the evap- 
oration process. As a result we have an increase in the 
radiation energy density, equivalent to a second phase of 
intense reheating, whereas the energy density in the PBH 
component falls down sharply. (We mention that there 
have been suggestions in the literature according to which 
the PBHs do not evaporate completely but leave behind 
Planck-sized remnants p^pTj . If this were the case, the 
final reheating would be reduced.) By direct inspection 
we confirm that there is no phase with a constant ratio 
between both energies densities (see the discussion pre- 
ceding and following Eq. (|52|) 



IV. THE EFFECTIVE ONE-COMPONENT 
MODEL 

In the present section we try to find an effective one- 
component description for the system of black holes and 
radiation. Since our considerations are based on the as- 
sumption that at some time to the temperatures T {BH) 
and T (2) coincide, it seems tempting to introduce an equi- 
librium temperature for the system as a whole, following 
the lines of [Q. The problem here is different, how- 
ever, compared with a mixture of two conventional flu- 
ids. Once the evaporation has started, an approximate 
equilibrium is no longer maintained, since there are no 
interactions like collisional events between the compo- 
nents. With their shrinking the black holes become hot- 
ter and hotter, while the radiation temperature decreases 
with the expansion of the universe, apart from phases 
of reheating. Therefore, the concept of an equilibrium 
temperature of the system as a whole intuitively might 
appear rather formal, except perhaps in the initial stage 
of the process. With the help of this quantity, however, 
a unified description of the cosmic matter will be possi- 
ble which comprises both conventional fluid aspects and 
properties characterizing the BH component. 

In a general two-fluid model close to equilibrium the 
Gibbs equation is 



Tds = d- + pd- - (p - p ) d— 

n n v 1 ' { ' n 



(57) 



where s is the entropy per particle. In the case of interest 
here the temperature T is the equilibrium temperature of 
the system of a conventional fluid and a PBH component. 
The temperatures T„, and T,„. of the previous sections 
do not appear as variables in the present effective one- 
temperature description. An (approximate) equilibrium 
for the entire system is usually assumed to be established 
through the interactions between the subsystems on the 



right-hand sides of Eqs. (||) and (Q). This reasoning 
does not hold in the case in which one of the compo- 
nents are evaporating black holes. Here, it is the circum- 
stance that radiation is evaporated at the BH temper- 
ature which may constitute a situation which is similar 
to a thermal equilibrium between two fluids. We assume 
that analogously to relations (|f^) the cosmic medium as 
a whole is characterized by equations of state 



p=p(n,n m ,T) , p = p(n,n {1) ,T) 



(58) 



If the first component is a black hole component, the first 
of the equation of state ( p8| ) specifies to 



P = P2 {n 2 ,T) 



(59) 



because of Eq. (]22 

If in the expressions ( |18[ ) for the entropies per parti- 
cle the temperatures T (1) and T (2) are identified among 
themselves and with T, and p (T) = p (1) (T) + p (T) as 
well as p (T) = p n ^ (T)+p (T) are used, the description 
based on Eq. ( |57j ) is consistent with the one relying on 
Eq. (|) for ns (T) = n (1)S(1) (T) + n (2)S(2) (T). 

The equilibrium temperature T is defined by f38| , ]39] | 

Pw ( n d) ' T (D ) + Pm ( n (2) ' T (2) ) = P ( n > n d) ' T ) ■ ( 60 ) 

As was shown by one of the authors |2^| , nonvanishing 
source terms in the particle number balances (^) 
(deviations from "detailed balance") give rise to a new 
type of effective bulk pressure for the cosmic medium as a 
whole ("reactive" bulk pressure), i.e. to enlarged entropy 
production. Here we try to apply this concept to the case 
in which the process responsible for the deviations from 
"detailed balance" is PBH evaporation. 

In order to find an explicit description we write the 
energy-momentum tensor of the system as a whole ten- 
tatively as §3] 



T' 



pu u 



( P + n) h l 



(61) 



i.e., we try to map essential features of the evaporation 
process on an effective bulk pressure n which we deter- 
mine below by consistency requirements. The relations 
= imply the energy balance 



p + e(p + p + n) = 



(62) 



for the effective one-temperature description. 

^From the Gibbs-equation (^) one finds for the change 
in the entropy per particle 

ns = 11 1 

T T 



n 



T 



[ r (1 ,- r (2,] • (63) 



The expression 
comes 



for the entropy production density be- 
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S° = nsT + ns 



(64) 



Introducing here Eq. (|63|) and the effective one- 
component chemical potential 



P+P rr 

fi = Ts , 



(65) 



we find 



sr. = -§n-fr 



"■(l) n (2) ( P (l) P(2) 

T 



[r (1) -r (2) ] . (66) 



Identifying again the first fluid with the black hole com- 
ponent, using Eqs. (|l7|), ( |2^ ) and (|3^) as well as the 
decomposition 

1 

n M = n (Bf f)Ai (B „) + n (2)M (2 ) = g'W) ' ( 67 ) 
(recall that p (2) = 0) one obtains 

s«=-|n-^r CBi0 . (es) 

The last expression for the entropy production density 
has to be consistent with Eq. ( ff"o| ) of the two-component 
description. This requirement provides us with the fol- 
lowing expression for the effective viscous pressure, 



9 



1 



m 



(BH) 



o P(B»|^(Bii| 



(BH) 
1 1 

f ~ f~ 

(2) 



1 



1 



T T 

^(BH) ^(2) 



(69) 



The last relation shows that the evaporation process may 
be described in terms of an effective viscous pressure II. 
In other words, PBH evaporation, if regarded as a spe- 
cific case of interfluid reactions, gives rise to a viscous 
pressure of the cosmic medium as a whole. The viscous 



pressure vanishes for T 



(BH) 



/to 



r 



o. 



i.e., no net change in the fluid particle numbers ("detailed 
balance"). This is also the limiting case of noninteracting 
fluids. We emphasize again that the quantity II in Eq. 
( |69| ) is exclusively a consequence of the evaporation pro- 
cess. "Conventional" bulk viscous fluid pressures, e.g., 
due to scattering of radiation particles by the BHs, have 
been neglected here for simplicity. 

By the field equations for a viscous cosmic medium 
with the energy-momentum tensor ( |6l] ) the quantity II 
influences the cosmic dynamics. Restricting ourselves 
again to the homogeneous, isotropic and spatially flat 
case we have 



87rp , 



-4tt (p + p + IT) 



(70) 



While physically the present one-component picture is 
equivalent to the two-component description of the previ- 
ous section, the mapping of certain features of the evapo- 
ration process onto an effective bulk pressure II may pro- 
vide a more transparent understanding of how the PBH 
evaporation modifies the cosmological dynamics. 

In Eq. ( |69| ) one has to deal with three generally dif- 
ferent temperatures: The temperatures T (BH) and XL. 
and the equilibrium temperature T of the system as a 
whole, defined by Eq. ( |60|). The dynamics of T {BH) and 
T . is given by the laws (|32|) and (|34|), respectively. But 
we have still to find a corresponding relationship for T. 
Such a law may be obtained via similar steps that led us 
to Eq. ([jjj). There exists, however, the following compli- 
cation. Because of the additional dependence of p on 
one has now three partial derivatives of p: (dp/dT) n n j , 
(dp/dn) T , and (dp/dn^j nT - The requirement that 
s is a state function now leads to ( j|0|, Q) 



dp 
dn 



ID 

n 



Tdp_ 
' ndf 



d_ 

df 



fa -n ) 



(71) 



generalizing Eq. (|il|), and the additional relation 



-M(D -M (2) ~ T §r -M (2) ) ■ (72) 



dnn 



(i) 



Using the Gibbs-Duhem relations 

d P {A) = V) s w dT ( A) + n w d P( A) ( 73 ) 
for T (A) = T together with Eq. Q, one finds 

(74) 



(A) dT 


_ P 
n 


A) 
(A) 




Consequently, the relations 1 


0) and © 


may 


as 








dp P + P T dp 






Pm 


dn n n dT 


n 




n <2). 


and 








dp _ P w 


_P 1 v 






"(2, ' 





(75) 



(76) 



respectively, with p {1) = p {1) (T) and p {2) = p {2) (T), since 
we are within the one-temperature description. 

Differentiating the second of the relations (Q) and ap- 
plying Eqs. @, ©, ©, ©, and ©, we obtain 



dp 
df 



f = - (e - r) t 



dp 
df 



[en + T( P + p)} 



[r w -r w ] . (77) 
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Using here p — p (BH} + p {2) and the second relation of 
Eqs. (|l7|), we find the evolution law for the equilibrium 
temperature of a system of a conventional fluid and a 
"fluid" of primordial single-mass black holes, 



dTT 



-r(e-r 



P{BH) 



j££) _ en 



in 



T 



dp 
df 



P 



(BH) 

r„ 



(78) 



The temperature law (|7gj comprises both the law ( p4D f or 
the fluid temperature and the BH temperature law 
This unifying feature may be considered as a justification 
of the equilibrium temperature concept for the cosmic 
medium. For p {BH) = 0, Eq. ( |7^ ) reduces to to the 
expression ( p4| ) for radiation while for p (2) = p {2) = II = 
0, equivalent to p — P {BH) and T = T {BH) , it coincides 



with the behavior (|32|), equivalent to Eq. (|28|). 

Inserting here the expressions ( [45|) at t = to, ( [46| ) and 
( p7| ) and taking into account II (to) = for T (BH) (to) = 
T (2) (t ) = T (t ) [cf. Eq. ©], Eq. @ at t = t reduces 
to the perfect fluid temperature law 



(79) 



This result again proves the consistency of the concept of 
an equilibrium temperature T of the system as a whole. 



V. ENTROPY PRODUCTION AND VISCOUS 
PRESSURE AT THE BEGINNING OF THE 
EVAPORATION PHASE 



The temperature laws (|32|), ( |34! ) and (|7£]) allow us to 
calculate the entropy production ( |68| ) and the viscous 
pressure (69) at the beginning of the black hole evapo- 



ration explicitly. For the temperatures we have in linear 
order: 



r(t)«T(to)-(t-*o)er^ 

op 



t (2) (t)«T w (t )-(t-t )(e-r (2) )T 

and 

T lBH) (to)-(t-to)T (Bff) -i 
Assuming 

T(t )=T (2) (to)=T iBH) (to) , 
the following relations hold at t — t : 



d P{2 , 



(80) 



(81) 



(82) 



(83) 



dp = P (2) - P(BH) Op = 1 P {2 ) 

dT T dp 3 p m - p (BH) 

P(BH) _ ry, P{BH) 



= T 



and 



dp/dT p - p 



dp dp {2) _ 1 p {BH) 



dp d P{2) 3 p - p 



(84) 



(85) 



Using the rates (45) at t = to, ( |46D and (}47|), we obtain 
t-t 1 



T : 



1 - 



3r 1-0 



-67 



T(t ) , 



t - 1 
3r 



(0t - P) 



T(t ) , 



(86) 



(87) 



and 



t_Jo 

3r 



T(to) 



Either from Eq. Q4Gp or Eqs. (|68|) and (|6J) we find for 
the entropy production density 



j. „ 



1 P(BH) 

2t T 



1 



t-tp 
9t 



(7 + 8/3 + 3H t) 



(89) 



The back reaction of the evaporation process on the 
cosmological dynamics is determined by the viscous pres- 
sure II which becomes in lowest order, 



n 



8/3 



1-/3 



1-/3 + 3H t 



ZHot 



t-tp 
18r 



(90) 



Except in the range 5/8 < /3 < 1, the quantity II has a 
negative sign because of relation (48). Obviously, /3 ss 1 
corresponding to p (2) (to) ~ P {BH) (to) is not a reasonable 
initial condition since II diverges. The linear approxi- 
mation breaks down for p {2) (to) ~ P {B h) (*o)- As follows 
from the first relation (|S4|), the quantity dp/dT changes 
its sign for p {2) (to) — P {BH) (^o)- A reasonable initial con- 
dition may be /3 = 1/2, i.e., p (2) (t ) = 2p (BH) (t ). In this 
case we find 



n 



i 



i 



6H t 



t-h 



'P(BH) 







(91) 



Since we restricted ourselves to the initial stage of the 
evaporation, t — t < r is valid. For (3 = 2 equivalent to 
P(bh) ( f o) = 2p (2) (t ) one obtains 



n 



n 

18 



3H r 



t-t 



-Pi 



((3 = 2) . (92) 



If the PBH component makes up almost all of the cosmic 
matter at t = to we have (3 3> 1 and 
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n 



1 - 



JL 

3H t 



t-t 



(BH) ' 



(/? » 1) . (93) 



In the opposite limit (3 <C 1 the viscous pressure is 



n 



5 

18 



1 



1 



3HnT 



t-tr 



'(BH) > 



(/? « 1) . (94) 



Although these explicit expressions for II are only valid 
in linear approximation, i.e. in the initial stage of the 
evaporation (we have J IX | < p {BH) in all cases, consistent 
with our initial assumption of small deviations from equi- 
librium), they provide an idea of the direction of the pro- 
cess. The role played by II in the held equations (70) may 
be thought of as a change of the "7-law" (7=1+ p/p) 
of the cosmic fluid according to 



7^7 



1 



n 



with 



n 



p 5-8(3 
p + p ~ 4 + /3 1-/3 



P + P 



1-P + 3H„t 



(95) 



t-t 
6r P{BH) 

(96) 



initially, i.e., 7 is effectively reduced. Via the field equa- 
tions, a negative viscous pressure tends to increase the 
expansion rate of the universe. 



VI. CONCLUSIONS 

Based on the remarkable fact that the BH temperature 
law r,_„. oc TYi~ naturally fits into the general formula 

(BH) (BH) J ° 

for the temperature of a fluid with variable particle num- 
ber, we set up a scheme in which the evaporation of PBHs 
is interpreted as a deviation from "detailed balance" in 
an interacting and reacting two-fluid system. A single- 
mass PBH configuration shares essential features with a 
pressureless, nonrelativistic fluid. We modelled the evap- 
oration process as a decay of such kind of "fluid" into 
a conventional relativistic fluid (radiation). We investi- 
gated the thermodynamics of this system, especially the 
entropy production and, with the help of the second law, 
derived general limits on the reheating of the radiation 
due to the evaporation. We found that intense reheating, 
i.e., an increase in the radiation energy density, may oc- 
cur both in the initial and final stages of the process. For 
our single-mass PBH ensemble there are no "equilibrium" 
solutions, i.e., solutions for which the ratio of the energy 
densities of the PBHs and the radiation remain constant, 
as obtained in ]2q| for a PBH configuration with a power- 
law mass spectrum. The impact of the PBH evaporation 
on the cosmological dynamics may be described in terms 
of an effective bulk pressure of the cosmic substratum as 
a whole which we evaluated explicity for the initial phase 
of the process. Its net effect is a higher expansion rate of 
the universe. 
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Figure 1.- Evolution of the energy density of the PBH 
fluid during the evaporation process. 



Figure 2.- Evolution of the energy density of the radiation 
fluid during the evaporation process. 
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